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higher-dimensional gauge theories. 
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I. INTRODUCTION 



The idea of unifying fundamental forces by introducing extra dimensions has attracted 
attention for many decades, and a theory realizing this idea is called Kaluza-Klein theory 

Recently, it has been observed by Arkani-Hamed, Dimopoulos and Dvali |2| that the 
existence of extra dimensions may play an important role to understand the hierarchical 
scales that exist between the weak and Planck scales. From a simple setting that only 
the graviton can propagate in the bulk corresponding to the extra dimensions while all the 
other fields of the Standard Model (SM) are located on a four- dimensional wall, they have 
concluded || that the length scale of the extra dimensions can be rather large ^ 10~ 2 cm, in 
contrast to previously suggested Kaluza-Klein theories in which the size of extra dimensions 
was of the order of the (four-dimensional) Planck length 10~ 33 cm or 1/Mqut ~ 10~ 30 cm, 
where Mgut is the unification scale in four-dimensional grand unified theories (GUTs). Their 
idea has been then followed and extended by several authors |3|,|]] to obtain more satisfying 
solutions of the hierarchy problem. Moreover, the above phenomenological proposal to 
confine fields on a lower-dimensional subspace fits well [§-|9| the D-branes || (extended 
objects attached by the end points of open strings) in string theories. 

If part of the SM fields can propagate in the bulk, and the size of the extra dimensions are 
large, the existence of such extra dimensions may be experimentally verified. There will be a 



number of phenomenological questions (see |T(| for instance) like " what are the experimental 
bounds on the size of the extra dimensions? JTTJ] " However, our concern in this paper is of 
a theoretical nature: Is the existence of a large extra dimension consistent with quantum 
theory? Our answer to this question will be "Yes" , provided that the compactification radius 
R is smaller than a certain value, the maximal radius Rm- It should be emphasized that 
the previous investigations [[I"2]-|14|| on non-abelian gauge theories in five dimensions on a 
lattice (which indicated that the theory have no continuum limit) were performed in the 
uncompactified case. These works fl"3|,[14|| were motivated to investigate whether or not the 
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non-trivial ultraviolet fixed point found in the e-expansion ]T5| is real. 



To be more specific, we consider pure SU (2) Yang-Mills theory in five dimensions where 
an extra dimension is compactified on a circle S with the radius of R. (It would be more 
"realistic" to compactly the fifth dimension on the orbifold S/Z2 so that the zero-modes 
contain only four- dimensional gauge fields and no scalar fields. We leave the case of S/Z 2 
to future work.) One may expect that the theory will carry the basic property of a four- 
dimensional gauge theory if the radius R is sufficiently small, while in the opposite limit of R 
the theory becomes more five dimensional. So there may be the maximal radius Rm below 
which the theory can possess a continuum limit with a non-vanishing string tension and 
can exist non-perturbatively. We will indeed find that our numerical simulations based on 
a compactified lattice gauge theory are supporting the correctness of this heuristic picture. 

The string tension is one of the most familiar physical quantities, which can give a 
physical scale to the lattice spacing. However, at a deconfining phase transition of first 
order, the string tension vanishes discontinuously, and we cannot use it for that purpose 
in this case. One of the crucial observations in this paper is that, if the fifth dimension is 
compactified, the first order phase transition changes its nature at a certain compactification 
radius. We will see this on anisotropic lattices by performing Monte-Carlo simulations with 
various compactification radii and by investigating the phase structure. The simulations 
also indicate that it could be possible to give a physical scale to the lattice spacing even in 
the deconfining phase if the theory is compactified, and this possibility will be studied more 
in detail. 

We will assume that the phase transition due to the compactification occurs at a certain 
value of R, the critical compactification radius Rc, and that the compactification radius is 
kept fixed at Rc along the critical line of the phase transition due to the compactification. 
That is, the critical compactification radius Rc is assumed to be a physical quantity. This 
assumption enables us to compute the lattice /5-function for a given R as a function of 
the lattice spacing 04 of the four- dimensional direction. In doing so, we can verify non- 
perturbatively the power-law running of the gauge coupling constant g, and find that the 



observed power law behavior fits well to the one- loop form suggested in Refs. |p|JT5|-p2| . The 
results for the lattice /5-function obtained from our Monte-Carlo simulations indicate the 
self-consistency of the assumption above. 

The results obtained for the lattice /3-function can also be used to make a further as- 
sumption on the physical scale in the deconfming phase and to investigate various scaling 
properties of the longitudinal Creutz ratio (defined in Eq. (|16|) ), making a discussion on the 
existence of continuum limits of the theory possible. We will be led to the interpretation 
that the theory may possess a continuum limit with a non- vanishing longitudinal string ten- 
sion if the compactification radius R is smaller than Rm ~ Rc/3, and that the non-trivial 
ultraviolet fixed point found in the e-expansion in the continuum theory may no longer be 
spurious. 

After we define our lattice action in Section II, we start to present the details of our 
calculations. In Section III we calculate the ratio of the lattice spacings £ = a^/a^ in terms 
of the parameters of the simulations (3 and 7, and then we discuss the phase structure in 
Section IV. In Section V we compute the lattice /3-function and then study on a continuum 
limit in Section VI, and the last Section is devoted for conclusion. 



II. THE ACTION 

In order to investigate the effects of a compactification in the 5-dimensional SU(2) gauge 
theory, it is crucial to employ an anisotropic lattice which has different lattice spacings, CI4 
and 05, in the four- dimensional directions and in the fifth direction, and is often used in 
the case of lattice gauge theories at finite temperature. We find that the effects of the 
compactification on an isotropic lattice can appear only for a small lattice size of the fifth 
direction (< 2) so that it is practically impossible to study the theory with different sizes 
of this direction. Another advantage is that, since we can vary 04 and independently, we 
can investigate the 04 dependence of physical quantities while keeping fixed. This enables 
us to study scaling properties in the compactifed theory for a given compactification radius 
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R. 

We denote the five-dimensional lattice coordinates by zm {M — 1, ... ,5) while the four- 
dimensional ones by x^ (fi = 1,...,4) and the fifth one by y. Link variable takes the 
form 

U M (x, y) = { U^(x, y) = U (x, y; x + a 4 fx, y) , U 5 (x, y) = U (x, y;x,y + a 5 ) } , (1) 

where U(z\ \ z 2 ) G SU(Nc) is the parallel transporter. The plaquette variables are 

Up 4 = U^(x, y) = U^x, y)U v {x + o 4 /t, y)C/J(ar + a 4 i>, y)Ul(x, y) 

Up 5 = U&(x, y) = U^x, y)U 5 (x + a 4 /2, y)Ul(x, y + a 5 )l/l(x, y) . (2) 

The Wilson action for pure SU(Nc) Yang-Mills theory in five dimensions is given by 

S = & E [ 1 - TT Re Tr U Pa ] + & £ [ 1 - i Re Tr U Ps ] , (3) 



/ J 4 



where J2p 4 — J2 Z i<^<^<4 an d Sp 5 = S 2 i<^<4- Periodic boundary conditions are imposed 
in all the directions. Q The coupling- and correlation-anisotropy parameters are defined as 

7 V/34 a 5 ' lj 

where 7 = £ is satisfied in the tree level. In the naive continue limit a 4 , a$ —>■ with the 
length of the fifth dimension fixed at 2irR, the action (|3]) becomes 

S = -E I(fj^) ^ + (^) ™& ] + 0(a 5 ) , (5) 



which goes to 



d 4 x ^ dy ~i TrF% IN , (6) 



if /? 4 = 2N c a 5 /g1 and /3 5 = 2N c a 2 Jg^a 5} where A M = # 5 ^M Ta > f mtv = d M A N - d N A 
i[ A M , A N ], and we have used 



M 



Another interesting case, i.e., orbifold boundary conditions which kill the scalar zero mode, can 
be archived by imposing U (x, y;x,y + 05) = U'(x, —y — 05; x, —y). 
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U„(x,y) = e i95aiA ^ , U 5 (x,y) = (*a*°*M*,v) . (7) 



On a lattice we mean a compactification if 



A4CI4 _ a A N 4 _ A4 

N 5 a 5 ~ 2irR ~ N$ > U 



is satisfied. Note that the gauge coupling constant g 5 has the dimension of y/al, and can be 
expressed as 

* 2 =2^'^V^- (») 
Later on we will use a dimensionless coupling constant g, 

g~> = (2nR)g? = ^ , (10) 

which is normalized for the four- dimensional Yang- Mills theory with the tower of the Kaluza- 
Klein excitations. At this point, Eq. fllOD is only a tree-level definition. 



III. f-y RELATION 

The parameters of the simulations are f3 and 7 for a given size of lattice, and the lattice 
spacings 04 and 05 are functions of these parameters. The introduction of an anisotropy 
into a lattice means that the regularization breaks 0(5) invariance of the continuum theory. 
To recover this symmetry we have to fine tune the anisotropy parameters 7 and £ that are 
defined in Eq. @. At the tree-level, it is £ = 7 as we have stated in the previous section. In 
higher orders the tree-level relation suffers from quantum corrections so that it can depend 
on (3 and 7, i.e., £ = £(7,/?). The basic idea to find the corrected relation, which has been 
intensively used in the study of QCD at finite temperature, is to use that symmetry. There 
are variants of the method, and we have decided to use a slightly modified method that is 



based on the matching of the Wilson loop ratio p3|-p5|. Let us briefly explain the method 
below. 
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We consider two kinds of Wilson loops W(zm, zn), the one W^x^jXu) within the four- 
dimensional subspace and the other one W(x^,y) that is extended into the fifth dimension, 
and calculate the ratios 

Since the Wilson loop is related to the static quark potential as 

W(zm, zn) ~ exp{-z M V(z N )} for z M -> oo , (12) 

we find that the rations ( |TT1) for large x and y become 

R(Xfj,,x u ) ~ exp{-a A V(x u )} , R(x^y) ~ exp{-a 4 V(y)} . (13) 

The 0(5) symmetry of the continuum theory requires then that 

R(x fl} x u ) = k R(x ll ,y) for x v = n u a A = y = n 5 a 5 , (14) 

where we have allowed the presence of the factor k. We measure the ratios for a given set 
of the lattice size, (3 and 7, and assume that they take the form 

R(x^, a^riy) ~ ki exp{— cr^n^} and R(x^, a 5 n 5 ) ~ /c 2 exp{— a^n^} , (15) 

and that they should become identical with each other, by symmetry, when n v a^ = 7150,5. 
From this consideration we obtain £ = 04/05 = o^jo^. Note that the ansatz ([X5|) has a 
meaning only in the confining region of the parameters, of course. 

In the practice, we fit the ansatz (|T^) for the data, and then scale by z (i.e. — > zn^) 
in such a way that -R(x M , za^n^) becomes closest to R{x lll 04^4), where we assume that k = 1 
on the r.h.side of Eq. ( |i~4|) 0. In the ideal case we would have z = o\joi = £. 



2 0n a lattice where one can obtain more data points, it is more convenient to use the method 
developed in Ref. 2B] for QCD, in which k is different from 1. In our case, due to the size of 



our lattice, we cannot obtain enough number of data points. In such case k = 1 is a reasonable 



assumption, as it has been discussed in Ref. |2J|. 



To restore the 0(5) symmetry in an efficient way, simulations are performed using the 
heat bath algorithm of the SU(2) lattice gauge theory. We employ a lattice of Nf x N 5 
shown in TABLE I, where ~ 7A4 is satisfied. We generate 5000 configurations, and 
Wilson loops are measured every 5 configurations. 

Fig. 1 shows £ versus (3 for various values of 7, and we see that £ is almost independent of 
(3. The data points for larger (3 are not plotted because they correspond to the deconfining 
region so that the ansatz ( |15|) has no meaning. The same data are plotted in Fig. 2 which 
shows the 7 dependence of £. The data are summarized in TABLE I. The central value of 
£ in the table is the average of the data points in Fig. 1 for a fixed 7. 

IV. PHASE STRUCTURE 

In this section we would like to investigate the phase structure of the five- dimensional 
theory defined by the action ([3]). It is known from the mean field analysis that higher- 
dimensional lattice gauge theories in more than four dimensions have a first order phase 
transition The studies of Monte-Carlo simulations [[12|JT3[1 also indicate that in the case 



of SU(2) gauge theory the first order transition occurs starting at D = 5. Our task is to 
extend these analyses to the compactified theory. To this end, we will be intensively using 
anisotropic lattices to take into account the compactification of the fifth dimension. 

A. Longitudinal Creutz ratio 

The string tension between two quarks that are separated in space is a typical physical 
quantity for the theory. What we know from experiments is that the string tension cr p hy S 
between two quarks that are separated in the four-dimensional subspace should be non- 
vanishing so that the potential between them is linearly increasing with the distance r. 



3 See for instance Ref. [26] and references therein 
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And the string tension is a good physical quantity to introduce a physical scale for other 
quantities obtained by a lattice calculation. If the underlying gauge theory is formulated in 
five dimensions, however, the feature of the linearly increasing potential is not automatically 
present, and in fact, the first order deconfining transition is found in Refs. |T^JI^] . 
We measure the Creutz ratio x(hj) defined as 

*<" ) = - to W,--i)vy ( i -!,,-) >■ (16) 

where W(i,j) is a rectangular Wilson loop with lengths of i and j. The Creutz ratio with 
large % and j becomes the lattice string tension cr lat in the case of the linearly increasing 
potential between two quarks. So, if a Creutz ratio with large % and j takes a non-zero 
value, the corresponding Wilson loop shows the area law which we regard as "confinement" . 
We consider the Wilson loops longitudinal to the four-dimensional subspaces, because we are 
interested in the confinement property in this subspace. We would like to demonstrate that 
the Creutz ratio behaves differently for different types of lattice. The results obtained from 
Monte-Carlo simulations on an isotropic lattice of size 8 5 (7 2 = 1.0) and on an anisotropic 
lattice of the same size (7 2 = 2.0 and 7 2 = 4.0) are shown in Fig. 3, where the vertical 
axis stands for the Creutz ratio, and the horizontal axis stands for f3 = vT^i/V We have 
generated 2500 configurations for each simulation point after thermalization, and Wilson 
loops are measured every 5 configurations for the calculation of the Creutz ratio. 

We see from Fig. 3 that the phase transition between the confining and deconfining 
phase exists around (3 = 1.64 in the case of the isotropic lattice (7 2 = 1.0) as it was found 
in Refs. [12,13] and around (3 ~ 1.73 and 1.77 in the cases of 7 2 = 2.0 and 4.0, respectively. 
We have performed the simulations starting with an ordered configuration with Um — 1 
(defined in Eq. ([!])) and with a disordered configuration, thereby obtaining clear hysteresis 
curves. The open symbols are the results of the ordered start and the filled symbols are 
those of the disordered start. Our results indicate that the transitions are of first order, in 
accord with the finding of Refs. for 7 2 = 1.0. 
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B. Transverse Polyakov loop 



In the uncompactified case, the Polyakov loop plays the role for an indicator of confine- 
ment. Here we consider loops which are transverse to the four-dimensional subspace and 
define the transverse Polyakov loop as 

^ = ^E^II^,y), (17) 

iV 4 x 1\ C y 

where z is a Z(Nq) phase factor (z Nc = 1) such that arg(L) G (—tt/Nc,tt/Nc)- In contrast 
to the longitudinal Creutz ratio flTED which we have discussed in the previous subsection, the 
transverse Polyakov loop ( |T7D has no direct physical meaning in four dimensions, because 
we do not identify the fifth direction with the temporal direction. We may say however that 
the quark currents running into the fifth direction are confined if the transverse Polyakov 
loop (L) vanishes. 

Fig. 4 shows the results of the transverse Polyakov loop on the 8 5 and 8 4 x 12 lattices for 
various values of 7, while, for comparison, the average of the plaquettes (lxl Wilson loop) 
for the same lattices is shown in Fig. 5. 2500 configurations have been used to measure the 
Polyakov loop and the plaquette for each point. As in the previous subsection, the open 
symbols are the results of the ordered start and the filled symbols are those of the disordered 
start. As expected, we obtain clear hysteresis curves, and so the transverse Polyakov loop 
and the average of the plaquettes also indicate that the phase transition is of first order. 



C. Compactification effects 

It may be worth pointing out that the compactified (D + l)-dimensional SU(Nc) lattice 
gauge theory belongs to the same universality class as the D-dimensional Z(Nq) spin model. 
The case of QCD at finite temperature T is a well-known example, where the temporal 
direction is compactified with the length T -1 . We expect the existence of a similar phase 
transition due to the compactification in our case, which is of second order, because the 
phase transition in the four- dimensional Z(2) spin model (Ising model) is of second order. 
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So, we repeat the measurements of the transverse Polyakov loop flT7|) and the average of 
plaquette for the compactified case. 

In order to take into account the compactification of the fifth-dimension, we use 
anisotropic lattices of size 8 4 x 4 and 8 4 x 6. The results for the transverse Polyakov loop 
with different 7 are shown in Fig. 6 and 7. (In Fig. 6 we have included the result on a 
12 4 x 4 lattice which shows that there are practically no finite size effects.) Noticing that 
the compactification radius R(= N^a^/2ii) becomes smaller for a given N$ as 7 becomes 
larger (see Fig. 2 and TABLE I), we observe that the nature of the phase transition changes 
due to the compactification. Namely, the interval of (3 in which two phases coexist becomes 
narrower as 7 increases, and there are no intervals for 7 2 2 for the 8 4 x 4 case and for 
7 2 4 for the 8 4 x 6 case, respectively. These phase transitions seem to be of second order. 
Observe also that the transition interval of (3 for r j 2 = 1.0 does not depend on N 5 , while, in 
contrast to this, the transition point f3c for the second order transition for a given 7 depends 
on N 5 . From these results, we conclude that the second order phase transition is caused by 
the compactification, and that the first order transition is not related to the compactifica- 
tion. In Figs. 8 and 9, we plot the average of the plaquettes for the 8 4 x 4 and 8 4 x 6 lattices. 
The results show that the transition becomes weak (like a cross over transition) starting at 
7 at which the first order transition of the transverse Polyakov loop turns to be of second 
order. (In Fig. 8 we have included the result on a 12 4 x 4 lattice to make it sure that finite 
size effects are negligible.) 

In Fig. 10 we show the qualitative nature of the phase structure in the /3 4 — (3 5 plane, 
which we have obtained from the result of this section. The "confining" and "deconfining" 
phases are separated by the critical lines of the first and second order phase transitions. The 
position of the critical line (bold line) of the first order phase transition does not depend 
on the lattice size, while that of the second order one (solid line) depends crucially on N$. 
Below the critical line in the $4 — /3s plane, the transverse Polyakov loop vanishes, and it 
is different from zero above the line. Note that this does not necessarily mean that the 
longitudinal Creutz ratio ( |16|) vanishes in the deconfining phase. The longitudinal Creutz 
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ratio fllBD corresponds to the "spatial string tension" in QCD at finite temperature, which 
is defined by the spatial Wilson loop, and indeed is non-vanishing even in the deconfining 
phase |27|]. Fig. 11 shows the longitudinal Creutz ratio versus (3 for the anisotropic lattice 
of size 8 4 x 4 with j 2 fixed at 4.0. The figure shows that the longitudinal Creutz ratio varies 
smoothly as (3 enters into the deconfining phase of the transverse Polyakov loop, indicating 
that it could be possible to give a physical scale to the lattice spacing even in that phase. 
Since indeed the spatial string tension is known to obey a scaling law at high temperature 



27|| , we may wonder whether some continuum limit in the present might also exist. The 
following subsections and sections are devoted to investigate this possibility from another 
point of view. 

In the case of QCD at finite temperature, the critical temperature Tq is a physical 
quantity. As in that case, it is well possible that the critical compactification radius Rq is 
a physical quantity, and that the lattice system on the different critical lines in the f3 — 7 
plane for different N 5 corresponds to the same physical system. As a first check, we estimate 
roughly the critical radius Rc for two critical lines of the second order phase transition at 
the end point. As mentioned (see also Fig. 16), at 7 ~ for N 5 = A and at 7 ~ 2.0 for 
N 5 = 6, the second order transition line merges in the first order transition line. The value 
of £ at the merging points, respectively, is 1.78 for 7 = a/2.0 and 2.78 for 7 = 2.0, where 
we have used the data in TABLE I. From the data on the Creutz ratio for the 8 5 lattice 
(Fig. 3), we find that the value of the longitudinal Creutz ratio at the transition points is 
approximately constant independent of 7, i.e. 

ciat = o-physal ~ 0.7 , (18) 

where we identify the longitudinal Creutz ratio xihj) with large i and j as the lattice string 
tension oi at . Using this, we find 



N 5 a 5c N 5 »- 12 



2tt 2tt£ c 



0.7 



Cphys 




0.30/^ ^ (ig) 

0.29/ v /cF^ 
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where £ c = a^/a 5c . These values are consistent with the assumption that the lattice system 
on the different critical lines corresponds to the same physical system. Eq. flI8|) also means 
that the value of a± at which the first order phase transition appears is approximately 
independent of 7, indicating that this value might have a sensible meaning. In the next 
section, we will do another check by using the lattice /3-function. 



V. THE LATTICE /3-FUNCTION 

We are interested in the physics in the four- dimensional subspace with a certain compact- 
ification radius. The anisotropic lattice we have used in the previous section is convenient 
for computations with different 04 while keeping the compactification radius constant. In 
this section we would like to compute the lattice /^-function in the four-dimensional subspace 
with the compactification radius R fixed at a certain value: 

Aat = , (20) 

dct4 

where g = g^/V^irR is the four-dimensional, dimensionless gauge coupling. We will cal- 
culate in the subsection B the /3-function for the compactification radius R at the critical 
compactification radius Rc using two lattices with different N 5 , where N 5 also corresponds 
to the number of Kaluza-Klein excitations. So, if the theory we investigate should be re- 
garded as a four-dimensional theory with only a few number of Kaluza-Klein excitations, the 
/3-function should depend explicitly on N$. On the other hand, if we obtain the same lattice 
/3-function for different N5, we are indeed dealing with a five-dimensional theory, and finite 
iVs or equivalently finite 05 effects may be regarded as negligibly small. First we would like 
to check this point. Another motivation is that we would like to examine non-perturbatively 
the celebrated power behavior of the running of the gauge couplings in higher dimensions, 
which we will use to give a physical scale in the deconfining phase of the transverse Polyakov 
loop and then to discuss the scaling behavior of the longitudinal Creutz ratio ([H]) in the 
next section. 
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Since the gauge coupling g and the lattice /^-function Aat are dimensionless, we may 
assume that the lattice spacings and 05 enter only in the combination £ = a 4 /a 5 . Fur- 
thermore, the perturbative analyses and also the discussion that follows below suggest that 
the correct variable is 

2nN 5 _ 2nN 5 a 5 _ (2n) 2 R 
s — - — — • (^1J 

£ 04 04 

This choice of the parameter has a non-trivial meaning: We may conclude that, if g really 
depends only on s, the continuum limit 05 — > with the compactification radius R fixed can 
be taken, and R can be regarded as a physical quantity in this sense. 

In the case of QCD at finite temperature, the critical temperature Tq is a universal 
quantity. The analogy for our case would be that the critical radius Rc is a universal 
quantity of the theory. So, the compactification radius would remain constant along the 
critical line in the (3 — 7 plane. However, there is a crucial difference compared with the 
case of QCD at finite temperature, because the critical lines in the present case merge into 
the region of the first order phase transition which is not related to the compactification. 
Therefore, this assumption is not reliable in the region in which the transition is of the first 
order. 

Keeping these circumstances in mind and defining the gauge coupling as 

9 - = arJfegf = ^ (22) 

on the critical line of the second order phase transition f], we can re-write Eq. ( p0|) as 

d In Rc 



Aat — Aat (1 , , I 

am 04 

A ^ r * 1 4£ [ 1 £ dp c d 1 ] (23) 



• 1 b 4 d 
with Aat = tt S T" 

iV 5 as 



Tc [s) 



£ d(3 c dj 



Pc d'j d£ 



N 5 (3 C 

where use have been made of Eqs. (||), @ and CpTOD . Here, we denote Aat for the (3 function 
with the assumption that the Rc is constant along the transition line. If there is no 04 
dependence of Rc, this assumption is correct so that Aat = Aat- 



This definition of the gauge coupling has the same form as the tree- level one (10) 



14 



Note that the critical lines in the (3 — 7 plane are different for different N$. In Eq. fl22|) 
we are implicitly assuming that g dose not depend on which critical line we use to calculate 
it. If we obtain the same gauge coupling from the different lines, it is a sign that the critical 
lattice systems for different N$ describe the same physical system. This will be checked in 
subsection B. 



A. Precise determination of the critical lines 

To compute the lattice j3 function /? lat using Eq. (p3| ), we need to know precisely the 
location of the critical points and its derivative with respect to 7 in the (3 — 7 plane. Let us 
therefore determine the critical lines in the /3 — 7 space next. To this end, we identify the 
transition point with the position of the peak of the susceptibility 

XL = iV 4 4 ((L 2 )-(L) 2 ), (24) 



where L is the transverse Polyakov loop defined in fll~Tp . We apply the histogram method 
|28| extended to an anisotropic lattice to evaluate the continuous parameter dependence 
of xl, as it was done in Ref. fl29f . To measure the Polyakov loop susceptibility, we take 



100,000 configurations. The results are plotted in Fig. 12 for iV 5 = 4 and in Fig. 13 for 
N 5 = 6. The large peak height at 7 2 = 2.0 for the 8 4 x 4 lattice and at 7 2 = 3.6 and 3.8 
for the 8 4 x 6 lattice (see Fig. 14) signals the first order transition which we have seen in 
the previous subsection. In Fig. 15, we see flip-flop in the history of the plaquette values, 
which is another sign for the first order phase transition. The transition point (3c and its 
derivative df3c/d'y for a given 7 are given in Table II. Here, the derivative of a transition 
point is calculated by fitting the continuous 7-dependence of (3c with the polynomial 

PcO)= E/n(7-7o) n , (25) 

ra=0 

where / n 's are fitting parameters, and d(3c/dry — fx at 7 = 70. The range of 7 and n max 
are chosen such that the results of the dPc/d-l are independent of the fitting range and the 
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fitting function. We adopt ±0.005 from the simulation point as the fitting range of 7 and 
the n max = 3 for the final results, respectively. The bin size of the jackknife error analysis is 
1000. 

The transition points in the f3 — 7 plane are shown in Fig. 16, where the circles (o) are 
the results for N 5 = 4 and the diamonds (o) are those for N 5 = 6, respectively. The short 
lines on these symbols denote the upper and lower bound of the slope of the transition curve. 
Two solid lines show the boundaries of the region in which two kind of phases coexist. Note 
that these boundary lines in Fig. 16 are obtained in the uncompactified theory. (Fig. 10 
is an illustration of Fig. 16 transformed into the ^4 — [3$ plane.) The interpolation curves 
are the dashed curves in Fig. 16, which are determined from the positions of ftc{l) an d its 
slopes. As we see from the figure, the critical lines bend strongly at (3 « 1.71 and 7 ~ 1.42 
for the 8 4 x 4 lattice, and (3 ~ 1.75 and 7 ~ 2.0 for the 8 4 x 6 lattice. The bending points 
are the merging points of two transition lines, the one for the phase transition characterized 
by the second order transition of the transverse Polyakov loop ([T7|) and the other one by 
the first order transition that is insensitive to N 5 . 



B. Calculation of /?i at 

Using the data given in Table I and II, we can express (3 function in terms of s, where s 
is given in Eq. (21). Then it is straightforward to compute (3\ at from Eq. (^3l). The results 



are shown in Fig. 17, where the o points are obtained on the critical line with N$ = 4 and 
the □ points are those with N$ = 6. As we see from Fig. 17, we obtain the same /^-function 
for two different N$ (or a^ c ). This implies that the lattice system on two different critical 
lines describes the same physical system, and finite N$ or equivalently finite effects may 
be regarded as negligibly small. In Fig. 18 we show g~ 2 defined in ( |22|) obtained from the 
data. This data indicate that g~ 2 depends only on the variable s, supporting our assumption 
that the critical compactification radius Rc is a physical quantity. Moreover, the Fig. 18 
suggests that g~ 2 {s) is almost a liner function. Its theoretical interpretation will be given in 
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next subsection. Note that the result above obtained for f3\ at does not verify the assumption 
that the compactification radius R is kept fixed at the critical value Rc along the line of 
the phase transition due to the compactification. To verify this assumption we need an 
analytical consideration as we will see in the next subsection. 

At this point we should emphasize that, in the region with small g 2 , the transition is of 
first order and is not related to the compactification. It implies that there is no reason to 
assume that the compactification radius is Rc near the first order phase transition. In Fig. 
18, the order of the transition turns to be of first order around g~ 2 = 0.95 for both cases of 
A^5 = 4 and 6. Therefore, the reliable region in which the compactification can be assumed 
to be Rc, is g~ 2 < 0.95. We, however, will assume in the next section, that the line of 
R = Rc exists, departing from the transition line around g~ 2 = 0.95 and entering into the 
deconfinement phase. How this line extends into the deconfinement phase cannot be found 
out within the framework of the Monte-Carlo simulations; we need analytical considerations 
as we will do in the next subsection. There we will discuss the theoretical interpretation of 
our data, and extrapolate the line of R = Rc into the region of a smaller g 2 . 



C. The e expansion, the power-law behavior and the ultraviolet fixed pont 

The power-law behavior of the gauge coupling is indeed suggested by its canonical di- 
mension, dim[g] = (4 — D)/2, where D is the number of the space-time dimensions. In 



the various explicit computations in perturbation theory |p|, |15|,|22| , this behavior has been 
directly seen. However, the explicit computations have been carried out basically within the 
frame work of perturbation theory, and so the result may not be trustful because the theory 
is perturbatively non-renormalizable 



5 The result of [ 18, 22 1 goes slightly beyond the perturbation theory because, though a number of 
non-trivial truncations to define an approximation scheme should be introduced, it is based on the 
exact Wilson renormalization group approach [pO]. 
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The simplest way to see the power law behavior in perturbation theory may be in the 
dimensional regularization scheme, as we do it briefly. Let go be the dimensionless gauge 
coupling in the pure SU(Nc) Yang-Mills theory in D = 4 + e dimensions. The /5-function 
is given by p"5||L3 | 



?D=A^f = eg 2 D + ^- 2 9 D + 0(e 2 ) with b = ~N C . (26) 
aA lbir z 6 

Now to mimic the dimensionless gauge coupling defined in the compactified theory (see 

(f22|)), we introduce 

~ g2 = (2ttRA) d ^ ' (27) 

whose /5-function becomes 

^ = A f = I^*A) D -V + ... (28) 

' 1 : (2ttRA) g* + ... for N c = 2 , D = 5 . (29) 



12vr 2 

Eq. ( pB"D suggests that there could exist a non-trivial ultraviolet fixed point for go, and as 

we have mentioned in Introduction, this possibility in the uncompactified theory was ruled 

out by the numerical studies of Refs. fli~2|-ri4}1 ■ Note that if the fixed point is real, then it 

means that the redefined coupling g behaves as an asymptotically free coupling, because 

- 2b (27rRA) D - 4 a , . 

g~ ^-— 2 p _ 4) ^oo as A-oo. (30) 

Translated into go-, we obtain 

-2 <T 2 2b 

9D = (2^)^^-^-4)16^ aS A ^°°' (31) 

which is consistent with the fixed point value obtained from Eq. ([26|). 

The form of the lattice /3-function in perturbation theory may be derived from the f3- 

function (p9"D, if we know the relation between A and a^. Since all the (four-dimensional) 

momenta in a lattice theory are restricted to the first Brillouin zone (— (it/ 0,4) < p M < 

(71-/(24)), the momentum cut-off is 1 7r / 04 1 . That is, A 2 = yj^ =1 (vr / 'a^) 2 = (27r/a,4) 2 , which 

implies that 



A=^. (32) 
a 4 



So, the suggested one-loop lattice /3-function is 

flffi = » 4 • (33) 

where f]we have used s = 2nN 5 /^, R = N^^tt and £ = a 4 /a 5 . 

D. The power law from the data 

Now we would like to proceed with our numerical analysis. Since the data in Fig. 18 
suggest that g~ 2 can be approximated by a linear function in the region we investigate, the 
one-loop form of the /3-function ([33]) is approximately correct. So we fit the function of s 
with the form 

for the data of g~ 2 at the second order transition. We find that the best values for N 5 = 4 
are C x = -0.208(8) and C 2 = 9.84(17) with X 2 /df = 1.4, and those for N 5 = 6 are 
C\ = —0.263(15) and C% = 10.71(21) with x 2 /df = 0.4, respectively. In Fig. 18 we compare 
the data for g~ 2 with g~ 2 for iVs = 6. As we can see from Fig. 18, the one-loop ansatz (|34"D 
fits well to the data, and moreover, the coefficient in front of s on the right-hand side of Eq. 
(|34]) is close to the one suggested in Eq. (|33|). Since the data with N 5 = 4 and 6 seem to 
lie on the same line, we also fit these data simultaneously. We obtain C\ = —0.224(6) and 
C*2 = 10.16(11) with x 2 /df = 1.7, which is a reasonable value, implying that the fitted gp 2, s 
for different A5 agree with each other. The fact that our data have a one-loop interpretation 



6 So far there exists no perturbative computation of /?i at in literature. Note also that the one-loop 
coefficient of depends not only on the regularization employed, but also on the definition of 



the gauge coupling. Our definition is given in Eq. (22) 
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indicate that the assumption that the compactification radius R is kept fixed at Rc along 
the line of the phase transition due to the compactification may be correct. 

Next, to discriminate the logarithmic behavior we would like to try to fit for the date on 
g~ 2 a function of the form 

gr 2 = B l + ^- 2 \ns, (35) 

and find that B x = -0.836(14) and B 2 = 100.9(1.1) using the data for N 5 = 4 and 6. 
This fit is not a good one because we obtain x 2 /df = 33. Moreover, the coefficient B 2 for 
the logarithmic function (|35|) cannot be explained within the frame work of perturbation 
theory. Namely, if the compactified theory on a lattice is simply a four- dimensional theory 
with Kaluza-Klein excitations of a finite number n, then the coefficient B 2 should be equal to 
(40/3)n. Since n could vary between 1 and N§ = 6, perturbation theory for this assumption 
would predict 

13 < B 2 < 80 , (36) 

which clearly disagrees with the value of B 2 obtained from fitting for the data. Since we have 
found that the one-loop form of the power law behavior describes the data well, the higher 
order contributions, especially those coming from non-renormalizable operators (remember 
the naive continuum theory is not renormalizable by power counting) must be suppressed, 
at least in the parameter region in our numerical simulations. 

It is the subject of the next section to investigate this possibility, where we will assume 
that the theoretical function (^) can be used to draw the lines of R = Rc even in the 
deconfining phase of the transverse Polyakov loop (|T7|). 



VI. TOWARD A CONTINUUM LIMIT 

In the weak coupling regime, which is the most important regime to investigate a con- 
tinuum limit, the way to use the transverse string tension as a physical quantity to give a 
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physical scale is not available, because the phase transition due to the compactification in 
that regime disappears. One of the central assumption in discussing the continuum limit in 
this paper is that the one-loop function (0) can be extended into the weak coupling regime. 
Equivalently, we assume that the /3-function for a given compactification radius R is given 
by 

dg 2 C 2 R 4 , , 

Put = ~a A1 — = — - g , (37) 

where Ci is given in Eq. (|34|). The assumption implies that we can draw lines of R = const, 
in the weak coupling regime. On these lines (3 becomes a function of 7, which can be obtained 
from TABLE I, (f22|) and fl34|) . So, given the lattice /3-function fl37|) we now know how to 
approach continuum limits. The object whose scaling property should be investigated is the 
longitudinal lattice string tension o"i at which we replace by the Creutz ratio fll6|), because, 
as we have seen in section IV, the longitudinal Creutz ratio can be non-vanishing even in 
the deconfining phase of the transverse Polyakov loop. In the following subsections we will 
investigate the scaling law of the longitudinal string tension 

O'lat = 0"phys a 4 = ^phys* 3 ^ 2 , (38) 

where we have assumed that <J p hy S should remain finite in the continuum limit. 



A. 05 — > limit 

We apply the scaling hypothesis ( |38|) to the 0,5 — > limit with £ fixed at a certain value. 
As we have stated, we assume that the one-loop function ([34]) can be extended into the weak 
coupling regime. If we move along the line of £=const., we change the compactification radius 
R. To express this more precisely, we first derive the scaling law for this case. To this end, 
let us consider the lines of R = const, for various A5 in the parameter space (/3,£), where it 
is implicitly assumed that the constants C\ and C2 in Eq. ([34]) are independent of A5 (the 
consistency of this assumption is checked for N 5 = 4 and 6 in section V): 
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D 1C 
(3 = --±£ + D 2 , D 1 = -4C l , D 2 = -l . (39) 

iV5 671 



Since a 5 = R c /(2irN 5 ) and Rc is assumed to be a physical quantity, we obtain from (|39| ) 

a 5 oc (iVs)- 1 oc (J3 - D 2 ) . (40) 

Inserting Eq. fl40|) into fl38|), we find that 

In (jjat = 2 In |/3 — D 2 \ + const. (41) 

Here, we use the value of D 2 obtained by fitting the data for N$ = 4 and 6 simultaneously. 
The result of this scaling behavior is shown in Fig. 19, where we have used the 8 4 x 4 lattice 
with 7 = 2.0 as for Fig. 11. The bold line corresponds to the theoretical line (f|l|). We see 



from Fig. 19 that the scaling law (il) is well satisfied for 1.6 ^ (3 ^ 1.8. Below ~ 1.5 we 
enter into the region of the strong coupling, and, above ~ 1.9, finite size effects due to small 
a 4 presumably start to become visible. So we may conclude that the data are consistent 
with the scaling law (f4"ip. This is an important result, and is indeed the only result which 
supports the correctness of the assumption that the lattice spacing has a physical scale even 
in the deconfining phase of the transverse Polyakov loop and of our way how to extend the 
lines of R = const, into that phase: This is an evidence for the existence of the fixed point 
suggested in the e expansion. 



B. a 4 /i? -> limit 

Next, we would like to investigate the scaling behavior of the longitudinal lattice string 
tension in the a^/R limit with R kept fixed. Since 2ttR = N5CI5, the lattice spacing 05 
is kept fixed in this limit for a given N 5 . Then, the string tension should obey the scaling 
law, 

0iat oc £ 2 or lncri a t = 21n£+ const. . (42) 

We compute on 8 4 x N$ lattices with N 5 = 2, 3, 4, 5, 6 and 8 the longitudinal Creutz ratio 
x(hj) along the theoretical line of R = const, on which iVs = 6 is critical. To determine 
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this line, we used the C% and C2 in Eq.(|34]) obtained from the data of N5 = 6. Note that for 
a given N 5 the compactification radius R is N 5 a 5 /(2n) = (N 5 /6)Rc- In Figs. 20, we plot 
lnx{i,j) as a function of ln£. If the slope of the lnx{i,j) is equal to 2, the scaling relation 
of Eq. ( f42"D is realized. In the N§ = 8 case, the results of the ordered start and disordered 
start are split, and the longitudinal Creutz ratio xihj) with large i and j of the ordered 
start fall drastically when we move from a large £ to a small £. This is in accord with our 
expectation, because the lattice system corresponds to the uncompactified. In the N 5 = 6 
case (the compactification radius R is equal to Re) the longitudinal Creutz ratios also start 
to fall down around ln£ w 1.1. Therefore, the lattice system above does not correspond to 
any four-dimensional theory, rather it describes a full five- dimensional theory. Keeping this 
in mind, we continue to consider the N 5 = 5,4,3 and 2 cases. 

The results are also shown in Fig. 20. As we see from these figures, the longitudinal 
Creutz ratios no longer fall drastically. Comparing the slope of these longitudinal Creutz 
ratios with the straight lines of the slope 2, we find that the longitudinal Creutz ratios for 
N 5 <; 3 decrease faster than £ 2 as £ decreases. If this continues to smaller £, i.e. smaller 0,4, 
we may conclude that that for N 5 ^ 3 the string tension cr p h ys decreases as a 4 decreases so 
that (Tphys vanishes in the a^/R — > limit. 

As we have observed, the slope of the Creutz ratios becomes milder as R decreases. This 
tendency of the milder-becoming slope with decreasing R is a real effect and not an effect 
of a finite N$, at least N 5 > 4 or equivalently R > (2/3)i?c, which we may conclude from 
the fact that our data show that the critical lattice systems with iV^ = 6 and 4 describe 
the same physical system. Some finite size effects may be present in the case of N 5 = 3 
and 2 in Fig. 20. Nevertheless, the tendency can be seen for these cases, too. It is this 
tendency of the milder-becoming slope that suggests the existence of a continuum theory 
with a non-vanishing string tension. If we assume that in the present case of setting the 
longitudinal Creutz ratio starts to scale according to the scaling law (|42| ) from N§ = 2 on, 
we obtain the maximal compactification radius 
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below which the compactified theory with a non-vanishing string tension could exist non- 
perturbatively. 

We would like to notice that, though the qualitative nature of the milder-becoming slope 
is real, the scaling behavior of the longitudinal Creutz ratio itself is sensitive to the choice 
of the extrapolation function fl34| ) that describes the lines of R = constant. It is therefore 
clear that for a more definite conclusion more refined analyses with a lager size of lattice are 
indispensable. 

C. "Simulated" N 5 -> oo limit 

To consider the iV 5 — > oo limit with R = const., we have to enlarge the size of our lattice. 
Instead of enlarging the size, however, we can simulate the limit with the data that we have 
already at hand. We would like to argue below that the second limiting process, a&jR — > 
with a 5 fixed at an a 5c (see Fig. 19) can be interpreted as an a 4 , a 5 — ► limiting process 
with R fixed, (as — > with R fixed is the same as N 5 — > oo with R fixed.) 

We have been assuming that the theoretical function given in Eq. fl39|) describes a set of 
the lines of R = const, in the (3 — £ plane for different N 5 . All lines so obtained are assumed 
to be physically equivalent: To each point on a line, there exists an equivalent point on each 
line. It follows then that all the points on a line described by Eq. (|39|) for a given N 5 can 
be transformed into a line that is parallel to the (3 axis in the (3 — £ plane. The mapping 
can be easily found, because the values of the gauge coupling g on the physically equivalent 
points should be the same. Since (3 does not change if the ratio £,/N 5 is fixed (see Eq. (|39"D ), 
the value of g does not change if we move along a line parallel to the £ axis (see Eq. (|2~2"D ). 
That is, to find a set of physically equivalent points we just have to move parallel to the £ 
axis. Therefore, moving along a line described by Eq. (^) for a given iVs can be assumed 
to be physically equivalent to moving along a line with £ = const, while changing (3 and 
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We can simulate enlarging N$ without changing N$. Since the compactification radius R is 
assumed to be a physical quantity, it remains unchanged during the transformation. 

Consequently, the scaling behavior of the Creutz ratios studied in Fig. 20 can be rein- 
terpreted as the scaling behavior along a line with £ and R kept fixed, where the scaling 
law appropriate for this limiting process is given in Eq. fl4"I|): The vertical axis ln£ in Fig. 
20 should be replaced by In |/3 — D 2 \ — 0.97 and the straight line should be understood as 
ln£ = 2 In \ (3 — D%\ + const., where we have used Eq. ([39]). We arrive at the same conclusion 
as in the a&/R — > case, which we do not repeat here again. But as we have stated there, 
the tendency of the milder-becoming slope with decreasing R is a real effect, at least for 
R > (2/3)Rc- This is so here, too, because our data show that the critical lattice systems 
with N 5 = 4 and 6 describe the same physical system so that the above mentioned trans- 
formation at least between the N 5 = 4 and 6 lines is trustful. The simulated iV 5 — > oo limit 
we have considered here should be regarded as a prediction of the real limit, at least for 
R > {2/3)R c . 

VII. SUMMARY AND CONCLUSION 

Our motivation in this paper has been to see, within the frame work of the lattice gauge 
theory, whether or not the non-trivial fixed point found in the e-expansion in the continuum 
theory of the pure SU(2) Yang-Mills theory in five dimensions is spurious in the case that 
the fifth dimension is compactified. We have used intensively anisotropic lattices to take 
into account the compactification. We have found that the compactification changes the 
nature of the phase transition: A second order phase transition, which does not exist in the 
uncompactified case, begins to occur, and turns to be of first order at a certain point. 

Under the assumption that the compactification radius R remains constant fixed at the 
critical value Rc along the critical lines of the phase transition due to the compactification, 
we have computed the lattice /3-function j3\ at , and found that /5i at as a function of s, obtained 
from the critical line of N 5 = 4 and 6, is the same (see Fig. 17). We have also found that 
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the gauge coupling on these critical lines is the same (see Fig. 18). From these observations 
we have concluded that the critical lattice system with — 4 and 6 describes the same 
physical system, and we are led to the assumption that this is the case for all N$. 

As we can see from Fig. 18, the power- law running of the gauge coupling (the solid 
line) is consistent with the data, which has a simple one-loop interpretation. This is the 
fact that supports the correctness of the assumption, at least for N$ = 4 and 6, that the 
compactification radius R remains constant fixed at the critical value Rc along the critical 
lines of the phase transition due to the copmpactification. 

At this point it is the natural thing to extend our findings into the deconfining phase of 
the transverse Polyakov loop: We have assumed that the lattice spacing has a physical scale 



even in the deconfining phase and the one-loop ansatz (|34|) can be used to draw the lines 



of R = const, in that regime. The investigation of the scaling law ([IT]) for the longitudinal 
Creutz ratio ([16]) shown in Fig. 19 supports the correctness of this assumption. At this 
stage, the existence of the non-trivial fixed point suggested in the e expansion might be 
evident. 

We have investigated the scaling behavior of the longitudinal Creutz ratio in the a^/R — > 
limit with R kept fixed, and found that the slope with which the Creutz ratios fall in the 
a^/R — > limit becomes milder as R decreases (see Fig. 20). From this tendency of the 
milder-becoming slope we are led to the interpretation that the compactified theory having 
a non-vanishing string tension could exist non-perturbatively if the compactification radius 
R is smaller than the maximal compactification radius Rm- Our estimate is: Rm ~ Rc/3 ~ 

It is clear that to make our interpretation more solid, we need not only refined and 
detailed numerical analyses but also analytical investigations. We hope that further studies 
will clarify the problems on the quantum realization of the old Kaluza-Klein idea. 
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TABLES 
TABLE I. £ — 7 relation 



7 2 




/3-range 


lattice size 


1.50 


1.438(57) 


1.51868 - 


1.66565 


8 4 x 16 


2.00 


1.784(50) 


1.55563 - 


1.69706 


8 4 x 16 


3.00 


2.340(40) 


1.59349 - 


1.73205 


8 4 x 16 


4.00 


2.779(34) 


1.60000 - 


1.75000 


8 4 x 16 


5.00 


3.161(39) 


1.65469 - 


1.74413 


8 4 x 16 


6.00 


3.490(33) 


1.61666 - 


1.76363 


8 4 x 20 


8.00 


4.062(39) 


1.62635 - 


1.76777 


8 4 x 24 


10.00 


4.617(35) 


1.50208 - 


1.73925 


8 4 x 24 


16.00 


5.923(51) 


1.50000 - 


1.70000 


8 4 x 32 
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TABLE II. Results for /3c and dj3c/^ by the histogram method. The simulations are per- 
formed at (/?4, fa). 



lattice 


7 2 


&) 


dc 


d/3 c /d 7 


8 4 x 4 


2.0 


(1.21250, 2.42500) 


1.71472(6) 


-0.0087(58) 




2.1 


(1.18350, 2.48535) 


1.71342(25) 


-0.067(16) 




2.5 


(1.07080, 2.67700) 


1.69060(36) 


-0.2219(82) 




3.0 


(0.95000, 2.85000) 


1.64702(31) 


-0.3337(67) 




4.0 


(0.77000, 3.08000) 


1.54018(48) 


-0.4049(79) 




6.0 


(0.55100, 3.30600) 


1.34560(59) 


-0.426(17) 




8.0 


(0.42500, 3.40000) 


1.19790(47) 


-0.368(11) 
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FIG. 1. (3 dependence of the anisotropy parameter £. 
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FIG. 2. Relation between £ and 7. 
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symbols are the results of the ordered start and filled symbols are those of the disordered start. 
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FIG. 4. Expectation values of the transverse Polyakov loop on an 8 5 lattice for 7 2 = 1.0, 2.0, 3.0 
and 4.0, and those on an 8 4 x 12 lattice for = 6.0. 
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FIG. 5. Expectation values of the plaquette on an 8 5 lattice for 7 2 = 1.0,2.0,3.0 and 4.0, and 
those on an 8 4 x 12 lattice for 7 2 = 6.0. 
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FIG. 6. Expectation values of the transverse Polyakov loop on an 8 4 x 4 lattice. The star 
symbols are the results on a 12 4 x 4 lattice. 
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FIG. 8. Expectation values of the plaquette on an 8 4 x 4 lattice. The star symbols are the 
results on a 12 4 x 4 lattice. 
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FIG. 10. Illustrations of the phase structure for the non-compactified case (left) and the 
compactified case (right), where a is the longitudinal Creutz ratio and L is the transverse Polyakov 
loop. 
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FIG. 12. [3 dependence of the Polyakov loop susceptibility obtained by the histogram method 
on an 8 4 x 4 lattice with 7 2 > 2.1. The circles denote the simulation point. 
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FIG. 13. (3 dependence of the Polyakov loop susceptibility obtained by the histogram method 
on an 8 4 x 6 lattice with 7 2 > 4.0. The circles denote the simulation point. 
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FIG. 14. Large peaks of the Polyakov loop susceptibility obtained by the histogram method 
at 7 2 = 2.0 on an 8 4 x 4 lattice, and 7 2 = 3.6 and 3.8 on an 8 4 x 6 lattice, respectively. 
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FIG. 15. Flip-flop in the history of the plaquette value at = 4.0 on an 8 4 x 4 lattice. 
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FIG. 16. Phase transition points for N$ = 4 (o) and N$ = 6 (O) in the (3 — 7 plane. Two solid 
lines denote the boundaries of the region in which two kind of phases coexist. Compare the figure 
with Fig. 10. 
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FIG. 17. jfliat as a function of s determined on the transition lines of N$ = 4 (o) and N$ = 6 
(□). The figures shows the physical equivalence between the critical lattice systems. 
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FIG. 19. Scaling behavior of the Creutz ratio along the j 2 = 4.0 line on an 8 4 x 4 lattice. The 
solid line is ln(x) = 21n|/3 - D 2 \ + 1.2 (see Eq.©), where D 2 = 2.156(23). 
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FIG. 20. Scaling behavior of the Creutz ratio for N$ = 2, 3, 4, 5, 6 and 8 measured on the line 
expected from the one-loop /3-function (|34|). The open symbols are the results of the ordered start 
and the filled symbols are those of the disordered start. 
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